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ABSTRACT 

Brightness distribution of Gamma-Ray Bursts (GRBs) is studied in detail under the 
assumption that GRB rate is related to cosmic star formation rate. The two populations 
of the long- and short-duration bursts in the 4B BATSE catalog are analyzed separately. 
Taking account of current uncertainties in the observational estimate of star formation rate 
(SFR), we have tried various models of the cosmic star formation history and we find that the 
SFR evolution in 2 = 0-1 is strongly constrained by the GRB distribution if the standard 
candle approximation is valid. The strong SFR evolution by a factor of ^ 15 from z = to 
1 inferred from UV observations is too steep to be consistent with the GRB distribution for 
any distance scale of GRBs. Some possibilities to reconcile this discrepancy are discussed, 
including the intrinsic luminosity dispersion of GRBs and/or modification of star formation 
history estimated by UV observations. We argue that SFR increase factor from z = to 1 
may be as low as about 4 if we choose different sets of cosmological parameters and/or take 
account of the evolution of metallicity and dust extinction in the UV data, and this would 
significantly remedy the discrepancy. 

Subject headings: binaries: close — stars: neutron — cosmology: observations — galaxies: 
evolution — gamma rays: bursts 



1. Introduction 

The brightness distribution of gamma-ray bursts (GRBs) observed by the Burst and Transient Source 
Experiment (BATSE) has been known to be significantly deficient in faint bursts compared with that ex- 
pected in the Euclidean space (Meegan et al. 1992), and this has been considered as one of the evidences 
of the cosmological origin of GRBs (Mao & Paczyhski 1992; Piran 1992; Dermer 1992). On the other hand, 
since most of GRB models are associated to death of massive stars whose lifetime is much shorter than the 
cosmological time scale, the brightness distribution of GRBs reflects not only the cosmological effects but 
also the cosmic star formation history (Totani 1997, hereafter T97; Sahu et al. 1997; Wijers et al. 1998). 
The cosmological origin of GRBs is now confirmed by the discovery of metal absorption lines in the optical 
counterpart of GRB 970508 (Metzger et al. 1997), and hence more detailed analyses of the GRB brightness 
distribution are required to investigate the origin of GRBs and the cosmic star formation history. 
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The most important consequence of the possible relation of GRB rate and star formation rate (SFR) is 

that the distance scale and absolute luminosity of GRBs become larger than those for no-evolution sources, 
because we need stronger cosmological effect in order to cancel out the rapid increase of the observed star 
formation rate from 2; = to 1 by a factor of more than 10 (Lilly et al. 1996). However, different authors give 
quantitatively different results: T97 pointed out that the scenario of binary neutron-star mergers (Blinnikov 
et al. 1984, referred to "the NS-NS model" , hereafter) results in a better fit than the case that GRB rate is 
simply proportional to SFR (referred to "the proportional model", hereafter), because the time delay during 
the spiral- in phase of binary neutron stars make the SFR evolution significantly flatter in z = 0-1. On 
the other hand, Wijers et al. (1998) concluded that the proportional model is in good agreement with the 
observed GRB brightness distribution when the redshift of the most distant GRBs is about 6. Petrosian 
& Lloyd (1997), however, claimed that neither the NS-NS nor the proportional model is consistent with 
the observationally determined SFR evolution. This paper investigates the origin of this discrepancy paying 
special attention to uncertainties in the SFR observations. We found that the GRB rate evolution in z = 
0-1 is crucially important in the GRB distribution analysis, and constraints on the cosmic star formation 
history in this redshift range are obtained, under the assumption of the standard candle approximation. 

It is well known that the duration of GRBs shows a bimodal distribution suggesting the existence of 
two populations of long ( > 2 sec) and short ( < 2 sec) GRBs (Kouveliotou et al. 1993). All previous 
papers including T97, which investigated the GRB brightness distribution in the context of the cosmic star 
formation history, however treated the GRBs as a single population. Because the energy spectrum of short 
GRBs is significantly harder than that of long GRBs (Kouveliotou et al. 1993) and the assumed spectrum 
affects the GRB brightness distribution analysis, we should analyze the GRB distribution separating the 
two populations. We define the long GRBs as those with Tqq > 2 sec, while the short GRBs with Tgo < 2 
sec, where Tgo is the interval over which 5 % to 95 % of the burst counts accumulate (Kouveliotou et al. 
1993), and analyze them separately in this paper. We also propose a new method to determine the average 
GRB spectrum in which the curvature of GRB spectrum compared to a pure power-law is taken into account 
based on the hardness-brightness correlation seen in the BATSE catalog. 

The paper is organized as follows: in §2, formulations used in this paper are described. We determine 
the average spectrum of the short and long GRBs and give the relations between redshift and peak flux of 
GRBs. In §3, the models of cosmic star formation history used in this paper are described. The results of 
flts to the observed GRB distribution of the 4B BATSE catalog (Paciesas et al. 1997) are given in §4. After 

discussing the result and its implications on the cosmic star formation history in §5, we conclude this paper 
in §6. Unless otherwise mentioned, we use the Einstein-de Sitter (EdS) universe with (/i, rioi^A) = (0.5, 1, 
0) as background cosmology. 



Expected number of GRBs in a flux interval between Pi and P2 {Pi < P2), where P is peak photon 
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Formulations 



flux in cm ^ sec ^, is given by 




(1) 



where dV/dz is the comoving volume element per unit redshift, ( the sky coverage efficiency 38 % for the 
BATSE experiment (Meegan et al. 1996)], Tobs the observation term, and Rqrb the comoving GRB rate 
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density. The factor (1 + z)^^ is included to account for the time dilation of the interval between detected 
bursts. We use the standard candle approximation throughout this paper and hence redshift and peak flux 
of GRB are related uniquely as z{P). The comoving volume element can be written as (c = 1): 



dV 

dz 



ATTdl 



dt 
dz 



(1 + ^) 

where t is the cosmic time and dt/dz in the Freedman universe is given as 

{ — ] = -Hoil + z)^{l + noz)il + z)2 - nAi2z + z2) 
and the luminosity distance d^ is given as 



dL = 



(l ~t~ z^dp^Qp , 



{Qo + f^A > 1) 

(r!o + nA = i) 

{Qo + f^A < 1) 



Here the proper distance dprop is given by 



(1 + ^) 



dz 



(2) 



(3) 



(4) 



(5) 



The most important fitting parameter in the GRB brightness distribution analysis is the distance scale 
of GRBs, and in this work we choose Zmax, the redshift of the most distant bursts corresponding to a threshold 
flux of the analysis, Pmin- The z-P relation is determined by the following equation: 

P _ l + Z cg|(Zniax) kp{z) 
^min 1 H~ ^max ^/^(^) ^p(^inax) 

where kp{z) is the /^-correction for the GRB photon spectrum: 

J Si de 

Here dLn/de is the differential photon number luminosity of GRBs at rest, and s the photon energy. The 
energy range in which the peak flux is measured, (£;,£„), is 50-300 keV for the BATSE experiment. The 
factor (1 + z)/(l + Zmax) in Eq. ^ is introduced because the standard definition of is valid for energy 
flux but P is photon number flux. The BATSE peak flux is given in three different time scales of 64, 256, 
and 1024 msec. We use the 1024 msec time scale for the long GRBs because the sensitivity is best in this 
scale, but 64 msec for the short GRBs because the average duration of the short GRBs is ~ 0.3 sec and 
we have to use shorter time scale than this for the real peak flux. We set the threshold values of Pmin = 
0.4 cm^^sec^^ (1024 msec) and 1.6 cm^^sec^^ (64 msec) for the long and short GRBs, respectively, above 
which the BATSE trigger efliciency is better than 90 % (Paciesas et al. 1997). There are 773 long GRBs 
and 192 short GRBs above these thresholds in the 4B catalog. 

Now let us determine the photon spectrum of GRBs. Because the spectrum is significantly different 
for the long- and short-duration bursts, we have to determine the spectra of these populations separately. 
Figure |l| shows a histogram of power-index of photon spectra (dLn/de oc e^") for the long GRBs (dashed 
line), short GRBs (dot-dashed line), and total (solid line). The values of a are estimated from the ratio of 
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thc flucncc in 50-100 keV and 100-300 keV in the 4B catalog. The average of a for long+short GRBs is 
1.38, which is consistent with a recent estimate of a = 1.1 ± 0.3 (Mallozzi, Pendleton, & Paciesas 1996). 
It is clear from this figure that the long-duration bursts have significantly softer spectra than short bursts; 
average a is 1.52 and 0.82 for the long and short GRBs, respectively. 

It is an easy option to analyze the GRB distribution with these average photon indices, but here we 
take a further step to include the effect of curvature of GRB spectra. Figure ^ shows a-P plots for the 
long and short GRBs. Both populations (especially the long GRBs) show a trend that weaker bursts have 
softer spectra, and this can be understood as follows. GRB spectra cannot be described by a pure power-law 
but exhibit some curvature, i.e., the spectrum becomes softer with increasing photon energy (e.g., Mallozzi, 
Pendleton, & Paciesas 1996). Because weaker bursts are more distant and hence more redshifted, the observed 
spectrum becomes softer for weaker GRBs on average. In the following we describe the determination of the 
GRB spectrum which takes account of this a-P correlation in a consistent way. First we assume a unique 
relation between a and P, which is obtained by a linear fit for the observed a-P relation, as shown in Fig. 
^ by the solid lines. We set a maximum peak flux Pmax, for which we adopt 60 cm^^ sec~^ both for the 
long and short GRBs, and denote the photon indices corresponding to Pmin and Pmax as apmin and apmax, 
respectively. We also assume a functional form of dLn/ds: 

log^^ = ^^(loge)^ +P log £ + const. (8) 

Although not physically meaningful, this form describes the four-channel burst spectra of the BATSE GRBs 
quite well (Mallozzi, Pendleton, & Paciesas 1996). The two parameters, A and B, are determined so that 
the average photon index in the rcstframe energy range of [(1 + z)ei, (1 + z)eu\ is consistent with the a-P 
relation: 

A(loge)p„^^ +P = -apmax (9) 

^(l0g£)p,„„ +B = -apmin , (10) 

where 

)ei] + log[(l + 

]} (11) 

(l0ge)p.„,„ ^ ^{l0g[(l-|-Zmax)ei] +l0g[(l + Zmax)£j} • (12) 

The above spectrum is determined if Zmax and Zminj which are the redshifts corresponding to Pmin and Pmaxj 
are fixed. The fitting parameter Zmax can be chosen arbitrary and the unknown parameter Zmin is determined 
by solving the following consistency equation: 



Pmax 1 ^" ^min ^/^(^max) ^p(-2'min) 

^/^(^min) ^p(-2max) 



(13) 



Now the GRB spectrum, which depends on the fitting parameter Zmaxj has been determined. We note that 
the above treatment of the GRB spectrum is the first to take account of the a-P correlation in a consistent 
way. In Figure ^, we show some examples of the GRB spectrum determined by the above method. The 
spectral curvature becomes larger with decreasing Zmax because the restframe spectrum has to become softer 
more rapidly with increasing photon energy. 

Now the relation between z and P is fixed, and the z-P plots are shown in Fig. ^ both for the long and 
short GRBs for some values of Zmax- Since the z-P relation has been fixed, we can calculate the cosmological 
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time dilation factor and redshift of some particular GRBs as functions of z,nax- Cosmological time dilation 
effect has been searched for various temporal structures, and a dilation factor of about 2 for the long duration 
bursts is widely accepted (see, e.g., Norris et al. 1994, 1995; Fenimore & Bloom 1995). Recently Che et al. 
(1997) suggested the existence of time dilation of about 2 also for the short GRBs. In the upper panel of 
Fig. ^, wc plot the time dilation factor, (1 + 2:wcak)/(l + -Zbright), as a function of Zmax- The brighter and 
weaker peak fluxes are taken as (fWeak, -Pbright) = (0.46, 9.5) and (1.6, 10.0) for the long and short GRBs, 
respectively. The values for the long GRBs correspond to the "dim+dimmest" and "bright" populations 
defined by Norris et al. (1994, see also Horack, Mallozzi, & Koshut 1996 for estimate of peak flux for these 
groups), and those for the short GRBs roughly correspond to the minimum and maximum of the brightness 
parameter defined by Che et al. (1997). The dilation factor of about 2 suggests Zmax ^ 2 both for the long 
and short GRBs. The lower panel of this figure shows the redshift of GRB970508, which is a long GRB with 
1024 msec peak flux of 0.85 cm^^sec^^ (when a = 1.5, see Kouveliotou et al. 1997), as a function of Zmax- 
The constraint of 0.835 < z < 2.3 for this GRB (Metzger et al. 1997) suggests 1.22 < z„iax < 3.43, although 
possible dispersion in GRB luminosity might relax this constraint. 

Although we have chosen the maximum redshift as the distance parameter of GRBs, it would be 
convenient to give relations between Zmax and absolute luminosity or total emitted energy of GRBs. If the 
gamma-ray emission is isotropic, the absolute peak luminosity in the restframe energy range of {ei,eu) is 
given as: 

/ \ 47r(i/^ (Zj^ax) -Pmin /-, -\ 

Vak = {£)t—F- UT1~7 T ' ^ ' 

where (e) is the average energy in restframe energy range of (ej, £„). The total energy emitted by GRBs in 
the same energy range can be estimated from /mim which is the fluence of GRBs at z = Zmax^ 

T-, 47r{i2v (^max) /min 

*°* ^ ft; — vTTT — ^ ' ^ ' 

where ke{z) is the /^-correction for the energy spectrum, 

Jsi ae 

We have estimated /min from fluence-P plots for the BATSE catalog. The peak luminosity and total energy 
emitted in the restframe energy range of 50-300 keV are shown in Fig. ^ as functions of Zmaxj as well as 
the values after bolometric correction (50-2000 keV). Since the spectrum becomes harder, the bolometric 
correction factor becomes larger with increasing z^ax- It should be noted that the bolometric correction is 
based on a simple extrapolation of the assumed functional form of dLn/de, but there might be cut-off energy 
in the GRB spectrum in reality. Therefore the estimate of bolometric luminosity or total energy has to be 
taken with care, especially for the short GRBs. 



3. Modeling the Cosmic GRB Rate History 

We can now calculate the GRB brightness distribution and compare it to the observation, provided 
that i?GRB is determined as function of z. In this paper we investigate two possibilities: the proportional 
model and the NS-NS model. Figure ^ shows observational estimations of comoving SFR density based on 
Ha or UV luminosity density (Gallego et al. 1995; Lilly et al. 1996; Connolly et al. 1997; Madau, Pozzetti, 
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& Dickinson 1998). In the Einstein-de Sitter universe with (fiojilA) = (1,0), the cosmic SFR at z = 1 is 
about 15 times higher than the local value, and the data at higher redshifts suggest a peak of cosmic SFR 
at z ~ 1-2. Let ^{z) be the ratio of SFR at z to the local SFR, i.e., ^(z) = 7?,(z)/i?*(0), where i?* is the 
comoving SFR density in the universe in [MQyr^^Mpc^'^]. In this paper, we parametrize the cosmic star 
formation history by two parameters of ^(1) and ^(4), and SFR evolution is determined by the natural cubic 
spline interpolation (e.g., Press et al. 1992) between SFRs at z = 0, 1, and 4. From the observational data, 
we adopt ^(1) = 13.6 and ^(4) = 3.43 as the standard observational star formation history (solid line in Fig. 

i- 

Although the UV luminosity density evolution in z = 0-1 measured by the Canada-France Redshift 
Survey (CFRS; Lilly et al. 1996) is very steep, another estimate of this evolution favors more modest 
evolution. Totani, Yoshii, & Sato (1997; hereafter TYS) pointed out that the present-day colors of galaxies 
are strongly correlated to recent star formation history in them and it is possible to estimate the SFR 
evolution at z < 1 from stellar population synthesis models of galaxies. TYS summed up the star formation 
history of galaxy evolution models for five types of galaxies, E/SO, Sab, Sbc, Scd, and Sdm, according to 
observed relative proportion of these types. Each of the five models is made to reproduce the spectral energy 
distribution in a wide range of wavelength at present and observed chemical properties (Arimoto & Yoshii 
1987; Arimoto, Yoshii, & Takahara 1992). The estimate of TYS for ^(1) is about 4 if the age of the universe 
is reasonable (10-15 Gyr). This value is significantly lower than the result of Lilly et al. (1996), but if 
one assumes a A-dominated universe with (rio,i^A) = (0.2, 0.8) the UV-estimated ^(1) becomes smaller by 
a factor of 1.8 and the discrepancy is a little remedied. Furthermore, Hammer et al. (1997) pointed out 
another effect which may be responsible for the very high ^(1) in the CFRS; star forming galaxies in the 
CFRS at z > 0.5 seem to have significantly low metallicities and hence they may be more transparent to 
UV light than local galaxies. This effect makes the evolution of UV luminosity density steeper than SFR 
evolution, and SFR evolution based on UV observations is overestimated when this effect is neglected. We 
will later discuss these points after we obtain the constraints on ^(1) from the GRB distribution. 

It is also widely discussed that the UV-estimated SFR at z > 1 is highly uncertain and may be 
significantly underestimated because of the dust extinction effect (Meurer et al. 1997; Pettini et al. 1997; 
Sawicki & Yee 1998; Heckman 1998; Cimatti et al. 1998). The ratio of far infrared luminosity to UV in 
metal-rich (~ Zq) local starburst galaxies is 10-100 (Heckman 1998) and this suggests that high redshift 
galaxies at z > 2, which also show signatures of starbursts, might be significantly affected by interstellar 
dust extinction by factor of ~ 10 (Meurer et al. 1997; Sawicki & Yee 1998) or in some cases ~ 500 (Cimatti 
et al. 1998). This possible upward correction for high redshift SFR data is interesting because the formation 
of spheroidal systems such as elliptical galaxies or bulges of spiral galaxies also predicts higher SFR at z > 2, 
if they were formed by starbursts at high redshifts. In order to give a quantitative estimate, let us calculate 
the total amount of star formation in elliptical galaxies. The total B-band luminosity density at present is 
1.2 X lO^hLse Mpc-3 (Lilly et al. 1996), where h = i?o/(100km s'^Mpc-i), with the fraction of elliptical 
galaxies of 0.28 (Totani, Yoshii, & Sato 1997). By using the mass-to-luminosity ratio {M/L)b ~ 16.2h for 
elliptical galaxies (Faber & Gallagher 1979), the total amount of stars which have been formed in elliptical 
galaxies becomes 7.6 xlO^/i^ [M© Mpc~^]. Here we have taken account of a recycling factor of 1.4 inferred 
from a galaxy evolution model of Arimoto & Yoshii (1987). The horizontal dotted lines in Fig. |^ are the 
average SFR expected in elliptical galaxies if they were formed before z = 2 or 3. This result suggests that 
the upward correction by a factor of ~ 10 is very likely from the viewpoint of spheroidal system formation. 
Considering the uncertainties discussed above, we will try various values of ^(1) and ^(4) in the following 
GRB analysis. We have plotted some examples of the model SFR evolution with different values of ^(1) 
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(=13.6 or 4) and ^(4) (=3.43 or 34.3) in Fig. |. 

If GRBs are associated to NS-NS mergers, we cannot use the simple approximation of i?GRB oc i?*, 
and the GRB rate evolution is written as: 



where t is the cosmic time and Pm{x) is the probability distribution of time delay from star formation 
to NS-NS mergers per unit mass of star formation, in units of [Mq^ yr"^]- The form of P„j can be well 
approximated as Pm{x) — A{x/tc)'' with tc ~ 0.02 Gyr and 7 ~ — 1 (T97). More detailed calculations 
of Pra by binary population synthesis suggest the uncertainty in 7 is probably about 0.2-0.3 (Lipunov et 
al. 1997; Portegies Zwart and Yungelson 1998). Because the aim of this paper is to investigate the cosmic 
star formation history from GRBs, we fix the parameters tc and 7 to 0.02 Gyr and —1, respectively. The 
normalization factor, A, can be inferred also from binary population synthesis calculations. Lipunov et al. 
(1995) estimated the present merger rate as 2 xlO^^ [yr^^], when 1O^^M0 stars have been formed constantly. 
Portegies Zwart & Yungelson (1998), by similar calculations, found the present rate as 3.4 xlO~^ yr~^, when 
the SFR is constant at 4 M0yr~^ and the present age of the Galaxy is 10 Gyr. Requiring these values to 
be consistent with the form of Pm, we find A = 1.6 xlO~^^ and 7 xlO^^"' [^^0^ yr~^] for the results of the 
two groups, respectively. We use A = 1 x lO"^'^ [-^^©^ yr""'^] for calculation of the absolute rate of NS-NS 
mergers and the uncertainty in A is probably about a factor of 2, as suggested from the difference of the two 
groups. This uncertainty affects the overall normalization of GRB rate and hence estimate of the beaming 
factor of GRBs, but does not affect the evolution of GRB rate which is the main topic of this paper. In 
Fig. 1^ we show some examples of calculation of Rqtib{z), generated from the four models of i?*(z) shown 
in Fig. ^ The GRB rate evolution becomes flatter than the star formation history when we take account 
of the time delay, and the importance of this effect is clear if one compares these results to those neglecting 
the delay (dotted lines). 



T97 used the Kolmogorov-Smirnov test for the comparison between the observation and the model cal- 
culation, which utilizes the cumulative flux distribution, A'^(> P). However, the differential flux distribution, 
dN/dP, reflects the GRB rate evolution better than A^(> P) and we use the analysis for the differential 
distribution in this paper. We set 9 bins in (Pmim^max) with a logarithmically constant interval both for 
the long and short bursts. The expected number of GRBs in each bin is calculated by using eq. (Q). 

Adopting the standard observational star formation history, i.e., ^(1) = 13.6 and ^(4) = 3.43, we have 
calculated the expected number of GRBs in each bin with nine values of Zmax (1^5 with an interval of 0.5), 
and the results are shown in the upper panels of Fig. || for the long bursts. (We have chosen the best-fit 
normalizations for the model calculations and hence only the shape of GRB flux distribution is tested.) In 
the upper-left panel the GRB rate is assumed to be proportional to SFR while the upper-right panel assumes 
the NS-NS merger scenario. These panels show that the inclination of the model curves is hardly changed by 
Zmax in peak flux larger than ^ 10 cm~^sec~^, while in small P range the expected event number becomes 
significantly smaller with increasing Zmax, because of the paucity of SFR at high z in the observational star 




(17) 



4. Comparison to the BATSE Data 



4.1. Brightness Distribution of GRBs 
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formation history. The standard observational SFR evolution does not agree with the GRB distribution, in 
either the proportional or NS-NS model. Because the inclination in the bright range of P is almost insensitive 
to ^max, we can set a stringent constraint on the SFR evolution at relatively low redshift (z < 1). Actually 
these results suggest that the SFR evolution inferred from the UV luminosity density at z < 1 (Lilly et al. 
1996) is too steep to be consistent with the GRB brightness distribution, both in the proportional model 
and the NS-NS model. One way to get a better fit is to decrease the value of ^(1). In the lower panels, a 
lower value of ^(1) = 4, as inferred from the galaxy evolution model (TYS), is used and then the discrepancy 
in the bright P range is significantly remedied. Figure |l^ is the same with Fig. ^ but for the short GRB 
population. Although the statistics is poor, the UV-estimated SFR evolution in z 0-1 is too steep also for 
the short GRBs if the GRB rate is proportional to SFR. 

The serious discrepancy between the UV-estimated SFR evolution in z = 0-1 and GRB distribution 
can also be seen in a SFR-z diagram. When GRB rate is proportional to SFR, we can convert the observed 
number of GRBs directly into cosmic SFR. Let A^obs(-pLj P2) be the observed number of GRBs in the interval 
oi Pi < P < P2. Then average SFR in the redshift interval of 2(^2) < z < z(Pi) can be written as: 

r^(Pi) dV dz ^ ' 

Jz{P2) dz (l + z) 

Figure ^ shows evolution of cosmic SFR calculated from the observed number of long GRBs as a function 
of z with various Zmax, where the binning in P is the same with that in Fig. ^. These curves are normalized 
to the observed local SFR at minimum redshifts. The SFR evolution inferred from GRBs becomes steeper 
with increasing Zmax but the evolution in z = 0-1 is much flatter than the UV-estimated evolution even 
when we take an extreme value of Zmax = 10. Therefore we conclude that the UV-estimated SFR evolution 
is inconsistent with the hypothesis that brightness distribution of GRBs traces the cosmic star formation 
history, if the standard candle approximation is valid. In order to get an acceptable fit, one must reduce the 
high SFR at z = 1, and this is consistent with the previous result of T97 who found that the star formation 
history based on the galaxy evolution model of TYS (i.e., ^(1) ~ 4) gives a good fit but the observational 
SFR evolution does not. This result is also consistent with that obtained by Petrosian & Lloyd (1997), but 
does not agree with that of Wijers et al. (1998). 



4.2. Constraints on Cosmic SFR Evolution 

In this section we give quantitative constraints on the cosmic SFR evolution, i.e., ^(1) and C(4), by 
detailed statistical comparison of the theoretical and observed brightness distribution. In the test, there 
are four parameters both in the proportional and NS-NS models: Zmax, ^(1), ^(4), and normalization of the 
GRB rate. We have searched the best- fit values which give the minimum chi-square (x^j^) in the two models 
of GRB rate evolution both for the long and short GRBs, and the results are listed in Table |l|. The searched 
ranges are 1-5 for z^ax and 1-100 for ^(1) and ^(4). The goodness-of-fit of the models can be assessed by 
X^in with 9 — 4 = 5 degrees of freedom. The significance level (S.L.), i.e., the probability that Xmin is larger 
than the observed value by statistical fluctuation is also given in the table, and large values of S.L. in the 
table show that the fitting is acceptable. 

Figure [l^ shows the allowed regions of ^(1) and ^(4) for the long and short GRBs, both for the 
proportional model and the NS-NS model. Here the allowed region of the ^(l)-^(4) plane is derived as 
follows. For each set of (^(1), ^(4)), we have searched the best-fit value of z^ax which gives a minimum value 
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of chi-squarc, which wc denote as XminKll), C(4)]- Then Ax^ = XminKll), ?(4)] — xLin is distributed as a 
chi-square distribution with 2 degrees of freedom (Press et al. 1992). Therefore the 95% C.L. region of (^(1), 
^(4)) is given by a contour of Ax^ = 5.99. It should be noted that the constraints on (f(l), ^(4)) obtained 
as above are conservative about the unknown value of Zmax- As discussed in the previous section, the most 
important result is the constraint on ^(1): for the long GRBs, ^(1) < 3.0 (proportional) and ^(1) < 3.8 
(NS-NS) at 95 % C.L., for any values of ^(4) and Zmax- The time delay in the NS-NS scenario makes the 
GRB rate evolution flatter in z = 0-1 and hence allows acceptable fits with larger ^(1) than the proportional 
model. The ^(1) is constrained also for the short GRBs as ^(1) < 4.5 and ^(1) < 10 (95 % C.L.), for the 
proportional and NS-NS models, respectively. If the SFR at z = 1 is more than 10 times higher than the 
present rate as suggested by the UV luminosity density, and if the standard candle approximation is a good 
approximation, the proportional model is completely rejected both for the long and short GRBs, and the NS- 
NS model is also inconsistent with the distribution of the long GRBs. In order to see the effect of changing 
the cosmological parameters, we show the same contour maps in Fig. |l3| only for the long GRBs in the 
two different cosmological models: open [(/i, JIq, J^a) = (0.6, 0.2, 0)] and lambda [{h, JIq, J^a) = (0.7, 0.2, 0.8)]. 
The constraint on ^(1) becomes weaker in the open universe, but the upper limit for ^(1) in the proportional 
model is still much smaller than 10. 

However, as discussed in §3, some independent estimates of ^(1) give lower values as low as ^(1) ~ 4, 
and there are acceptable models with such values of ^(1). In order to see the allowed region about z^ax in this 
case, let us consider a statistical test under the condition that a particular value of ^(1) is given. For a fixed 
value of ^(1) = 4, we have searched the best-fit values of (zmax,^(4)) and results are shown in Table |^ for 
the long GRBs with six model parameter sets (proportional or NS-NS, and the three cosmological models). 
The goodness-of-fit is acceptable except for the proportional model in the EdS and A-dominated universes, 
and allowed regions of (2:max,C(4)) for the acceptable models are shown in Fig. ^by thick contours. These 
contours are defined by the increment of x^ relative to Xmin i'^ Table || with two degrees of freedom. We 
have also repeated the same analysis but with ^(1) — 3, with which all the six models are acceptable (see 
Table and the results are given in Fig. |l^. The NS-NS model gives a better fit than the proportional 
model because of the flatter evolution due to the time delay. The maximum redshift is at least larger than 2 
and likely larger than 3, depending on ^(4). If the distance scale of GRBs is determined independently, the 
cosmic SFR at high redshift {z ^ 4) can be estimated. 

We have also calculated the best-fit normalization parameters, i.e., production rate of GRBs per unit 
mass of star formation (events/M©) for the proportional model and the beaming factor {Air/ Ail) for the 
NS-NS model, for these acceptable models. The results are shown by the thin contours in Figs, [ij and [T^ . 
The production rate of GRBs required for the proportional model is ~ 5 x 10~® events / Mq if the emission 
is isotropic, and the NS-NS model requires a beaming factor of ~ 200-400. 



4.3. Long and Short GRBs 

Because of the poor statistics of the short GRBs, we cannot derive strong constraints on them. However, 
it may be interesting to see whether there is a difference between the distance scale of the long and short 
GRBs, provided that the two populations are originated by the same physical event, i.e., the GRB rate 
evolution of the two populations as a function of z is the same. In Figure |l^ we plotted the reduced x^ of 
the fit as a function of Zmax for the long and short GRBs in the NS-NS model, using a fixed star formation 
history [(C(l);?(4)) — (4,34.3)]. This result suggests that Zmax of the long GRBs is larger than that of the 
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short GRBs. Combined with the result of Fig. this rcsuh confirms that the total energy emitted by the 
short GRBs in the restframe 50-300 keV is more than one order of magnitude smaller than that by the long 
GRBs if the beaming factor of these two populations is the same. 

It may be rather interesting to note that, although the total energy is different for the two populations 
by a factor of more than 10, the peak luminosity of the two populations is remarkably similar (see Fig. 
^). There may exist a physical mechanism which regulates the peak luminosity to a certain value in the 
shock dissipation of relativistic kinetic energy, and it may give a hint to understand the origin of the two 
populations. 



5. Discussion 

As shown in this paper, the star formation history in the universe estimated by UV observations is 
completely inconsistent with the brightness distribution of GRBs if the standard candle approximation is 
valid. The origin of the discrepancy is too steep evolution of SFR in the relatively low redshift range of z 
= 0-1. It is of course an option that GRBs are not related to death of massive stars, but considering that 
most of theoretical GRB models are associated to collapse of massive stars, possibilities of reconciling the 
discrepancy under the condition that GRB rate is related to SFR should be considered. 

Clearly one possibility is intrinsic luminosity dispersion of GRBs. Krumholz, Thorsett, & Harrison 
(1998) suggested that luminosity dispersion by a factor of more than 20 can reconcile the UV estimated SFR 
evolution and GRB distribution, assuming a power-law distribution of GRB luminosity. However, although 
some luminosity dispersion is likely as suggested from diversity of observational properties of GRBs, it is 
not clear that the luminosity dispersion is actually such a broad one. Krumholz et al. (1998) claimed the 
redshift estimations for GRB970508 {z = 0.835; Metzger et al. 1997) and GRB971214 (z = 3,42; Kulkarni 
et al. 1998) as an evidence of broad luminosity dispersion. However, as is well known, the absorption line 
at z = 0.835 gives only an lower limit for the redshift of GRB970508. The identification of the host galaxy 
for GRB971214 is also not perfect but there is a doubt that the 'host' galaxy is accidentally on the line of 
sight. The shape of distribution of GRB luminosity function is also unknown. The luminosity dispersion of 
GRBs is a possibility to reconcile the discrepancy but it requires further observational confirmation. 

On the other hand, it has recently been discussed that UV-estimated star formation history may be 
significantly biased due to neglecting the effect of interstellar dust extinction (Meurer et al. 1997; Pettini 
et al. 1997; Sawicki & Yee 1998; Heckman 1998; Cimatti et al. 1998), as discussed in §3. Therefore 
the possibility that the strong SFR evolution in z = 0-1 inferred from UV observations is overestimated 
should be considered. Hammer et al. (1997, hereafter H97) pointed out that CFRS galaxies at z ^ 1 show 
signatures of low metallicity and smaller dust extinction, which may cause an overestimate of SFR evolution 
in z = 0-1. H97 made a simple model to include this effect, in which the conversion factor between SFR 
and [Oil] luminosity depends on the restframe {U — V) color of galaxies and the dependence on the color 
is determined by SFR calibration of local galaxies. Then they found that the SFR evolution is ^(1) ~ 6.9 
from the CFRS data in the Einstcin-de Sitter universe, which is considerably lower than the pure evolution 
of 2800 A or [O II] comoving luminosity density (^(1) ^ 15). If we change the cosmological model into the 
open [(rio,f^A) = (0.2,0)] or the A-dominated universe [(^Oj^^a) = (0.2,0.8)], the value of UV-estimated 
^(1) is further reduced to ^(1) = 5.33 and 3.81, respectively. These values are now near the independent 
estimate (^(1) ~ 4) based on galaxy evolution models which reproduce the present-day properties of galaxies 
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(Totani, Yoshii, & Sato 1997). Therefore we conclude that more modest evolution of SFR in z = 0-1 is 
not unreasonable and this would significantly remedy the discrepancy between UV observations and GRB 
distribution. 



6. Summary and Conclusions 

In this paper we have presented a detailed study on the possible relation between the brightness 
distribution of gamma-ray bursts and the cosmic star formation history. The long and short GRBs in the 
4B BATSE catalog are analyzed separately. We proposed a new method to determine the average GRB 
spectrum in which the curvature of GRB spectra compared to a pure power-law is taken into account based 
on the observed a-P correlation. Various models of the cosmic star formation history are tried considering 
the present uncertainties in the observational estimate of SFRs, and implications of the GRB distribution 
on the star formation history and galaxy evolution were discussed. 

We have shown that the evolution of SFR in z = 0-1 is crucially important for the fit of GRB brightness 
distribution, and SFR evolution in this range is strongly constrained if the standard candle approximation 
is valid. The analysis on the long GRBs suggests that, in the Einstein-de Sitter universe with {h, ^Iq, I^a) = 
(0.5,1,0), the SFR increase factor from z = to 1 [^(1)] should be smaller than 3.0 (95 % C.L.) if the 
GRB rate is proportional to SFR, and than 3.8 if the GRBs are produced by binary neutron-star mergers. 
For short GRBs, ^(1) is constrained as ^(1) < 4.5 and ^(1) < 10 for the proportional and NS-NS models, 
respectively. These values are significantly smaller than the current estimate (^(1) ^ 15; Lilly et al. 1996) 
of SFR evolution based on the UV luminosity density, but marginally consistent with a theoretical estimate 
(^(1) ~ 4) by galaxy evolution models based on the local properties of galaxies (Totani, Yoshii, & Sato 
1997). These results are consistent with Totani (1997) in which the star formation history based on the 
galaxy evolution model gave a good fit while the observational history did not. 

We have discussed some possibilities to reconcile this apparent discrepancy between the UV-estimated 
SFR evolution and GRB distribution, under the condition that GRBs are related to death of massive stars. 
One possibility is intrinsic luminosity dispersion of GRBs, although neither the width of dispersion nor 
the shape of distribution is well known. We have also argued that the uncomfortably large ^(1) inferred 
from the UV observation may be an overestimation and the real value could be as low as about 4 if we 
choose different sets of cosmological parameters and/or take account of the evolution of metallicity and dust 
extinction. Therefore the UV observation itself can also be consistent with the BATSE data and the galaxy 
evolution model. 

For the case of low values of ^(1) ^ 4, we have obtained some constraints on the distance scale of 
GRBs and production rate of GRBs. We could not find any acceptable fit for the long GRBs with Zmax < 2. 
Therefore Zmax is at least larger than 2 and likely in the range of Zmax — 3-5, if the long GRBs are associated 
to death of massive stars or NS-NS mergers. The production rate of GRBs from star formation is ^ 5 x 10~^ 
[Mq^] if GRB rate is proportional to SFR, and the beaming factor required for the NS-NS merger scenario 
is about a few hundreds. 

The maximum redshift of the short GRBs seems smaller than that of the long GRBs if the GRB rate 
evolution is the same for the two populations. This confirms that the total energy emitted by the short 
bursts is smaller than that by the long bursts by more than one order of magnitudes (see Fig. if the 
beaming factor is the same. On the other hand, the peak luminosity of the two populations is remarkably 
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similar 2 x lO'^^ erg sec~^ in 50-300 keV). This might suggest that there is a physical mechanism which 
regulates the energy loss rate of shock heated matter in shocks generated by relativistic motion, in a wide 
range of total energy liberated as a relativistic fireball. 

The author would like to thank an anonymous referee for careful reading of this manuscript and useful 
comments. He has been supported by the Research Fellowships of the Japan Society for the Promotion 
of Science for Young Scientists, and the Grant-in-Aid for the Scientific Research Fund (No. 3730) of the 
Ministry of Education, Science, and Culture of Japan. 
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Table 1: Best-fit Parameters for 2;max, Ci^), and ^(4) 



ismology" 


Duration 


Model*- 






C(4) 


■ 

Amin 


S.L. ( 


EdS 


long 


prop. 


3.7 


1.7 


4.6 


2.34 


80 


EdS 


long 


NS-NS 


4.8 


1.6 


13.8 


0.85 


97 


EdS 


short 


prop. 
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5.48 
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short 


NS-NS 
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12.6 
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long 


prop. 
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1.7 


2.1 


2.00 
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long 


NS-NS 


4.7 


2.5 


18.2 


1.22 


94 


A 


long 


prop. 


2.-5 


1.0 


1.0 


2.05 


84 


A 


long 


XS-XS 


4,8 


1.4 


7,2 


I.IG 


95 



"Cosmological parameters used are (/i,Qo,Oa) = (0.5, 1, 0), (0.6, 0.2, 0), and (0.7, 0.2, 0.8) for EdS, Open, 

and A, respectively. 

''The proportional model in which GRB rate is proportional to SFR and the NS-NS model in which GRBs 
are produced by neutron star mergers. 
'^The degree of freedom is 9 — 4 = 5. 



Table 2: Best-fit Paraiiioters for "max ((f)- (Only for tlie long GRBs and ^(1) is fixed.) 
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NS-NS 
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2.3 


89 


4 


NS-NS 
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9.91 
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3 
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13 


3 
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25.1 


4.1 


67 


3 


NS-NS 


Open 


4.6 


21.9 


1.4 


97 


3 


NS-NS 


A 


4.5 


17.4 


5.6 


47 



"The proportional model in which GRB rate is proportional to SFR and the NS-NS model in which GRBs 
are produced by neutron star mergers. 

''Cosmological parameters used are {h,no,^A) = (0.5, 1, 0), (0.6, 0.2, 0), and (0.7, 0.2, 0.8) for EdS, Open, 
and A, respectively. 
'^The degree of freedom is 9 — 3 = 6. 
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Fig. 1. — The histogram of effective photon index of GRB spectra (a) in the 4B BATSE 
catalog (Paciesas et al. 1997). The photon index is calculated from the fluence ratio of the 
50-100 keV and 100-300 keV ranges. The dashed hne is for the long GRBs, the dot-dashed 
line for short GRBs, and the solid line for the total. 
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Fig. 2. — Correlation between the photon index (a) and peak flux (1024 msec for the long 
GRBs and 64 msec for the short GRBs), in the 4B BATSE catalog (Paciesas et al. 1997). 
The solid hnes are least-square fits. 
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Fig. 3. — The photon spectrum of GRBs used in this paper. In our treatment, the spectrum 
is determined when the fitting parameter, Zmax is fixed (see text). The curvature reflects 
the a-P correlation shown in Fig. ^ The upper three curves correspond to the long GRBs 
and the lower to the short GRBs. The solid, dashed, and dot-dashed hues are for the cases 
of -Zmax = 1, 3, and 5, respectively. The dotted lines are simple power- law spectra with the 
average photon indices of the 4B catalog [a = 1.52 and 0.82 for the long and short GRBs, 
respectively) . 
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Fig. 4. — The redshift vs. peak flux relation for the long and short GRBs, for some values 
of -2max which is the redshift corresponding to P1024 = 0.4 and Pq4 = 1.6 cm~^sec~^ for the 
long and short GRBs. The values are indicated in the figure. 
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Fig. 5. — Time dilation factor and redshift of GRB 970508, as functions of -Zmax- For the 
definition of time dilation factor, see text. In the upper panel, the solid line is for the long 
GRBs and the dashed line for the short GRBs. In the lower panel, the three lines are all for 
the long GRBs with different cosmological parameters: {h, Qq, Qa) — (0.5, 1, 0), (0.6, 0.2, 0), 
and (0.7, 0.2, 0.8) for the solid, dashed, and dot-dashed lines, respectively. The constraint 
of 0.835 < z < 2.3 which is set by the absorption lines (Metzger et al. 1997) is also shown. 
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Fig. 6. — Peak luminosity and total energy emitted by GRBs, as functions of the fitting 
parameter, -Zmax- The thick fines are for the long GRBs, in the energy range of 50-300 keV 
(solid) and 50-2000 keV (dashed) at the rest frame, while the thin fines are for the short 
GRBs. Isotropic emission is assumed in this figure. 



- 21 - 




Fig. 7. — Cosmic star formation history as a function of redshift. Data points are direct 
observational estimates by Gallego et al. (1995, open triangle), Lilly et al. (1996, filled 
circles), Connolly et al. (1997, open circles), and Madau et al. (1997, filled square). Curves 
are the model of cosmic star formation history used in this paper, with the parameters 

(^(1),^(4)) = (13.6, 3.43), (13.6, 34.3), (4, 3.43), and (4, 34.3) for the solid, dashed, short- 
dot-dashed, and long-dot-dashcd lines, respectively (sec text in detail). The two horizontal 
dotted lines arc average star formation rates expected if the stars in present-day elliptical 
galaxies were formed before 2; = 2 or 3. 
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Fig. 8. — The merger rate evolution of binary neutron stars as a function of redshift, cal- 
culated from the cosmic SFR evolution shown in Fig. |^. The four panels correspond to 
the different four models of cosmic SFR evolution with values of (^(1),^(4)) indicated in 
the figure. The line marking corresponds to that of Fig. |^. The dotted lines are the SFR 
evolution before convolved with the time delay, normalized at 2 = for comparison. The 
time delay between star formation and merger makes the rate evolution flatter m z = 0-1. 
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Fig. 9. — Differential distribution of peak ffux (1024 msec) of the long GRBs. The data 
points are the BATSE data and the error bars in peak; flux indicate the width of bins. Nine 
model curves are shown for each panel, corresponding to different values of ^max = 1^5 with 
an interval of 0.5. The number of GRBs at small peak flux decreases with increasing z^aa.^ 
because of the low SFR at high redshifts. The left panels are for the proportional model, 
while the right panels for the NS-NS model. In the upper panels we use = 13.6 while 
in the lower panels ^(1) = 4. In all panels ^(4) = 3.43. 
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Fig. 10. — The same as Fig. ^ but for the short GRBs with 64 msec peak flux. 
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Fig. 11. — The star formation history in the universe deduced from the brightness distribu- 
tion of the long GRBs (crosses combined with sohd hues). The GRB rate is assumed to be 
proportional to SFR. Four different values of 2;max are used, as indicated in the figure. The 
data points of filled and open circles, and filled squares are the observational SFR estimates 
as shown in Fig. ^ 



-26- 



100 




Fig. 12. — Constraints on the cosmic star formation history from GRB brightness distribu- 
tion. The allowed regions of ^(1) and ,^(4) are shown for the long and short GRBs (upper 
and lower panels, respectively) and for the proportional model and the NS-NS model (left 
and right panels, respectively). The sohd line shows the 95 % C.L. constraints, while the 
dotted lines for 68 % C.L. The Einstein-de Sitter universe with h — 0.5 is assumed. 
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Fig. 13. — The same as Fig. [T^, but only for the long GRBs with different cosmo logical 
parameters: {h,QQ,Q\) = (0.6, 0.2, 0) and (0.7, 0.2, 0.8) for the upper and lower panels, 
respectively. 
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Fig. 14. — Allowed regions in -2max-'C(4) plane are shown, assuming ,^(1) = 4. The left panels 
are for the proportional model and the right for the NS-NS model, with the three different 
sets of cosmological parameters for each, as shown in the figure. The thick-solid line is 95 
% C.L. constraint and the thick-dotted line for 68 % C.L. In the upper-left and lower-left 
panels, there is no allowed region. The thin solid lines are contours of the normahzation 
factor: production rate of GRBs per unit mass of star formation [GRB M^^] in the left 
panels and beaming factor (47r/AQ) in the right panels. 
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Fig. 15. — The same as Fig. |1^, but ^(1) = 3. 
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Fig. 16. — The reduced of the fit to the BATSE data as a function of Zmax, with a fixed 
cosmic SFR evolution [(^(1), ^(4)) = (4, 34.3)] in the NS-NS model. The solid line is for the 
long GRBs and the dashed hne for the short GRBs. The Einstein-de Sitter universe with 
h — 0.5 is assumed. 



